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Abstract 

We study in a quantitative way CP-violating effects in neutrino oscillation experiments in 
the light of current and future data. Different scenarios with three and four neutrinos are 
worked out in detail including matter effects in long baseline experiments and it is shown 
that in some cases CP-violating effects could affect the analysis of a possible measurement. 
In particular in the three neutrino case we find that the effects can be larger than expected, 
at least in long-baseline — > v e . Moreover, measuring these effects could give useful 
information on the solar oscillation frequency. In four neutrino scenarios large effects are 
possible both in the — *> v T and — > v e channels of long-baseline experiments, whereas 
short-baseline experiments are affected only marginally. 
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1 Introduction 



Understanding the mechanism responsible for the patterns and values of fermion masses, 
mixings and CP-phases is a very important goal of particle physics. The recent advances 
in the field of neutrino physics, especially by experiments measuring neutrino oscillations, 
are in this context extremely valuable. There is now overwhelming evidence for neutrino 
masses and it is possible to extract from data very interesting patterns for masses, mixings 
and as we will see CP-phases. Neutrino masses require extra ingredients, which constitute 
a small extension of the Standard Model and could be viewed as a complication, but it 
seems likely that this new information is an important second lever arm for the fermion 
mass problem. The reason is that in many extensions of the Standard Model, like for 
example in GUT theories, the apparent similarities between quarks and leptons lead to 
connections between lepton and quark masses via the same symmetry breaking Higgses 
and/or related Yukawa couplings. The new results from the neutrino sector can therefore 
often be combined with information from the quark sector allowing thus a better test of 
proposed mass structures (like e.g. in the form of so-called textures). Neutrino masses 
imply mixings in the lepton sector which include also CP-violating phases. These CP- 
phases are however not only interesting due to their appearance in mass textures, but they 
can be responsible for very important physical effects. Mechanisms have for example been 
proposed, where CP-violation in the neutrino sector is the essential ingredient in explaining 
the baryon asymmetry of the universe. CP-violation in the neutrino sector is thus by itself 
an important issue. 

Besides these questions of general theoretical kind the presence of CP- violation can also have 
quite significant impact on how experiments should be performed and analyzed. We will 
discuss the possibility of having CP-violation in various experiments from a qualitative and 
quantitative point of view and will point out that the effects can be larger than expected 
from qualitative arguments Jl|. In long baseline — > v e experiments we will find for 
example with three (four) neutrinos and maximal CP-phase CP-asymmetries up to 40% 
(60%), while short baseline experiments could see in some corners of parameter space of v^- 
v T oscillations effects up to 10% with four neutrinos, so that a modest improvement of the 
baseline would be enough to see sizeable effects. CP-violation can thus affect in a significant 
way (or even spoil) the measured probabilities, such that the planning of experiments, 
data analysis and theoretical interpretation should take them into account. Furthermore 
measuring CP-violating effects would provide information on the frequency giving rise to 
the solar oscillation, at least in the framework with — > v e transition and three neutrinos. 
In this respect CP-violation offers an interesting possibility to investigate in long-baseline 
experiments parameters typically accessible through solar neutrino experiments. This is 
because the solar frequency suppression of CP-violation is only linear and because a 6*13 
angle suppression affect the CP-conserving amplitude more than the CP-violating one. 

The minimally extended Standard Model with three heavy right-handed neutrinos allows 
usually only for two independent oscillation frequencies, while experiments claim three dif- 
ferent values. Therefore even complex three neutrino scenarios, where more than one pair of 
neutrinos oscillates, can not accommodate all three results simultaneously. Once this possi- 
bility is excluded one can either extend the theoretical framework such that all data can be 
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accommodated or one must discuss different scenarios excluding some of the data. We will 
consider both possibilities. Since the LSND evidence is considered most controversial we 
will exclude their result in an analysis with three neutrinos. The analysis of four neutrino 
scenarios is done with all available experimental information including LSND. In this case 
a third squared mass difference Attend m uch larger than the solar and atmospheric ones 
explains the LSND oscillation signal and can generate small CP-violation effects even in 
short-baseline experiments sensitive to this squared mass difference. 

Our paper is organized as follows. First we introduce our notation and framework. Next we 
give a general discussion of potential CP- violating effects in neutrino oscillation experiments 
and present the existing experimental data that we use in our analysis. Then follows an 
extensive analysis of CP- violating effects in scenarios involving three or four neutrino species. 
In the end all results are discussed in comparison and conclusions for planned neutrino 
oscillation experiments are drawn. 

2 CP-violation in neutrino oscillation experiments 
2.1 Masses, mixings and notation 

The three neutrino flavour eigenstates u ei , i 6 {1,2,3} form in the Standard Model (SM) 
with their respective charged lepton partners doublets under SU(2)l and without loss of 
generality we can define the neutrino flavour eigenstates in a basis where the charged lepton 
mass matrix is already diagonal. Neutrino oscillations require non-degenerate (Dirac or 
Majorana) neutrino masses, which are however forbidden in the SM, since any mass term 
or Yukawa coupling which generates neutrino masses would violate the gauge symmetry. 
Neutrino oscillations require therefore SM extensions with new fields and/or interactions. 
A new SU(2)l Higgs triplet field, for example, with hypercharge Y = 2 and non-vanishing 
vacuum expectation value can generate a 3 x 3 Majorana mass matrix Ml for the neutrino 
fields u ei via Yukawa couplings[|. Another possibility is the existence of further neutral and 
sterile?] neutrino states u e ,, j G {4, . . . ,un} resulting in the most general mass term in the 
Lagrangian 

£mass = (m N )ij U e M ej + h.C. , (1) 

with the extended x neutrino mass matrix 

(\M L M D \ ,„. 
m " = { MJ \M R ) ■ < 2 > 

In addition to Mi this symmetric mass matrix contains the 3 x [n^ — 3) dimensional 
off-diagonal entries Mp which can be generated by Dirac-like Yukawa couplings and the 

1 Note that for phcnomcnological reasons the VEV of such a triplet should be rather small compared to the 
electro-weak scale since it would contribute already at tree level to custodial SU(2) breaking. Furthermore 
the masses of the so far unobserved single and double charged Higgses should be very large. 

2 I.e. SU(2)l singlets. The existence of further non-singlet SU(2)l representations containing a neutral 
fermion is strongly disfavoured. Such representations would require further fermions to satisfy the anomaly 
conditions. Existing mass bounds for new fermions (e.g. generations) would then lead to big unobserved 
radiative corrections in the so-called S and T variables. 
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electro-weak VEV. Note that the remaining (un — 3) x (n^ — 3) dimensional sub-mass matrix 
Mr can be made of explicit mass terms with arbitrary values, since none of the symmetries 
of the Lagrangian is in this way violated. The elements of Mr are thus uncorrelated with 
the electro- weak sector and a natural range for these mass terms is the scale where the 
new fields become members of multiplets in some extended framework^. This implies that 
some (but not necessarily all) eigenvalues of Mr are heavy enough to decouple after the 
diagonalization of the matrix in eq. (Q), thus leaving only n > 3 light states involved in the 
low energy physicsQ. 

Independently of the physics giving rise to them, we will discuss in this paper the case of 
three or four light neutrino degrees of freedom, i.e. n = 3,4. These mass eigenstates are 
assumed to be mixtures of the original three flavour eigenstates which couple to SU (2) i and 
possibly with other sterile states u ei , i e {4, . . . , n}. In the limit where all heavy neutrino 
degrees of freedom are decoupled we can thus write down a n x n neutrino mass matrix 
which leads upon diagonalization to n = 3, 4 mass eigenstates and a n x n CKM matrix 
Uckm- Neutrino oscillation is in this picture to a very good approximation the oscillation 
of n neutrino degrees of freedom, where only one of the first three flavour eigenstates can 
be produced and detected (via its coupling to W's). These flavour eigenstate can be writ- 
ten with the help of Uckm as a superposition of the n mass eigenstates. The transition 
probability becomes then 

(3) 

where n x n matrix multiplication is implied, D = Diag(e~' lElt , . . . , e~ lEnt ) and where the 
indices i,j correspond to the respective flavour eigenstate with i,j e {1,2,3}. 

The n x n CKM matrix Uckm contains as usual a number of global unphysical phases 
which can be absorbed into the fermion fields. Because of the potential Majorana nature 
of the neutrinos more physical phases survive, however, compared to the quark case with 
pure Dirac masses. Altogether there are up to n(n — l)/2 physical mixing angles and up 
to n{n — l)/2 physical phases. For n = 3 (n = 4) we have thus when all masses are 
non-degenerate and not purely Dirac-like three (six) mixing angles and three (six) physical 
phases. In the n = 3 case we use the standard parameterization 

Uckm = U(6 12 , # 23 , 6 13 , 6) Diag(e ia \ e ia \ 1) , (4) 

where the so-called CKM-like phase 5 and the extra Majorana phases a\ and a 2 are ex- 

3 For tin = 6 the three new neutrino fields can for example be placed very economically together with 
evi, i S {4, 5, 6} in doublets of SU(2)r in left-right symmetric models. Contributions to Mr arise then via 
left-right symmetric Yukawa interactions resulting in masses of the order of the left-right breaking scale. 
Alternatively the new neutrino fields can be fitted into representations of some GUT group and natural 
entries in Mr would in this case be at the GUT scale. 

4 Another reason why some neutrino states might decouple is that there are two degenerate eigenstates 
one of which can be made decoupled with a rotation in their subspace. This would be the case for example 
if only the non-diagonal blocks in eq. (0) were non- vanishing (pure Dirac case). 
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plicitely shown. With the help of = cos(%) and s^- = sin(#y) we can chose as usual: 

C12C13 S12C13 
U(9 12 , 9 23 , 9 13 , 5) = \ -s 12 c 23 - c 12 s 23 s 13 e lS c 12 c 23 - s 12 s 23 s 13 e tS s 23 c 13 |. (5) 

S12S23 — Ci 2 C 23 S\ 3 e l5 —Ci 2 S 23 — Sl2C23Sl3e l<5 

For n > 3 we will not use any particular parameterization, but we can still factorize 
Uckm = U ■ Diag(e* ai , . . . , e"*"- 1 , 1). Since only U C kmDU^ km is involved in oscillation 
experiments, and since D is diagonal and commutes with the extra diagonal Majorana 
phases Diag(e* ai , . . . , e* 0n_1 , 1), it is always possible to eliminate those extra phases. Thus 
only {n 2 — 3n + 2)/2 phases (one for n — 3, three for n — 4) show up in oscillation exper- 
iments just like in the quark case. Due to this similarity we call these remaining phases 
"CKM-like". 

2.2 CP-violation 

The oscillation probability for a neutrino produced in a flavour eigenstate v e . to be detected 
as a u e . after having traveled a distance L with a ultra-relativistic energy E is 

P{y e% ^>v e .) = P C p{v e . ^v ej ) + P&p(v ei -> v ej ) , (6a) 

P{p e% ^Ve,) = Pcp{y ei -> i/ e .) - Pfi^(z/ ei -> i/ ej .) , (6b) 

where 

Pcp (^ -> = Sij - 4 ^ MJX 1 ) sin2 A kh , (7a) 

Pop K-^ e ,)=4£ Im( J fc e f ) sin A fc/l cos A fch , (7b) 

fc>h 

with J fc e f = U^Ut^uM^ Jkh = (JX 1 )* = {JhkY and A kh = Am^L/AE. Notice 
also the relations Ylk^kh, = ^ij\Uih\ 2 , Yli^kh^ = $kh\Ujh\ 2 - One can easily see that only 
n — 1 out of the n(n — l)/2 frequencies |Amj^J, k > h, are independent. From CPT 
invariance follows in general P[y ei v ej ) = P{v ej — > ^ e J and in particular -P(f ei — > z/ e J = 
P(z/ 6i — >■ z/ e J. Therefore one can see that CP-violating effects can not occur in disappearance 
experiments. 

Finding sizeable CP-violating effects in oscillations is however not easy since they are af- 
fected by suppressions similar to the quark case. First \Pqp\ < Pcp, since Pqp + Pqp 
and Pqp — Pqp have both to be positive. Moreover, the CP-violating contribution to the 
oscillation probability is suppressed, because CKM-like CP-violation is not possible with 
only two neutrinos (or only two non-degenerate neutrinos). As a consequence CP- violating 
effects need (at least) three mixing angles between non-degenerate neutrinos such that the 
squared mass difference corresponds to wavelenghts neither too large compared with the 
distance travelled (otherwise the oscillation does not have enough time to develop) nor too 
short (otherwise the effect is washed out). These requirements are made explicit by the 
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three neutrino version of eqs. (||[7|). Using cr^- = J2k £ ijk an d ^Jkh 1 = ~ a ij a khJcv it is in 
fact in this case 

Pc p (y ei -* u ej ) = Sij - 4 Re sin 2 A 12 - 4 Re J% sin 2 A 23 - 4 Re j£ sin 2 A 31 , (8a) 
Poppet -> ^) = ScTijJ sin A 12 sin A 23 sin A 31 

= 8(Tij J(sin 2 A 23 sin A i2 cos A i2 + sin 2 A i2 sin A 23 cos A 23 ) , (8b) 

with 

8 J = cos sin(2fl 13 ) sin(2#i 2 ) sin(2# 23 ) sin 5 . (9) 

From eqs. (^|b,|9|) we see explicitely that even in case of maximal CP- violation (i.e. | sin5| = 
1) small angles and small sinAi 2 suppress the effect. On the other hand it is important 
to notice that these two suppressions are only linear so that in cases in which it is safe 
to neglect sin 2 A 12 effects in the CP-conserving part of the probability, still CP-violating 
effects proportional to sinAi 2 can be relevant. Experiments sensitive to CP-conserving 
effects with amplitudes suppressed by the sin 2 of a small angle, can thus have a chance to 
see a CP-violation effect only proportional to the sin of that angle. We will study this in a 
quantitative way in the following. 

2.3 Asymmetries 

As a measure of CP-violation, we will consider suitable asymmetries between CP-conjugated 
transitions. Besides having obvious physical meaning these asymmetries will show to what 
extent the analysis of a possible signal in a single channel v ei —>■ u e . (or v ei — > p e .) performed 
without taking into account CP- violation could be spoiled by CP- violation effects. Neutrinos 
travel in some experiments through matter such that the two conjugate channels have to 
be distinguished in the vicinity of the corresponding MSW region. Where appropriate we 
will use therefore the suffix "m" to express that the transition probabilities in eqs. @ are 
changed by the presence of matter. With this in mind we define the total asymmetry 

tot = (P m (^^^))-(P m (^^^)) 

where the average symbol () in eq. (|T0|) accounts for the averaging in energy and length 
present in every real experiment and is particularly important in case both the channels, 
and therefore the asymmetry itself, are measured. In this case 

<P m (^ - v ej )) = f d(L/E)P™(v ei - v ej ){L/E)f{L/E) , (11a) 
<P m (P es - u ej )) = J d{L/E)P™{v ei - u e .)(L/E)f(L/E) , (lib) 

where the weight functions /, / include the initial spectra, the cross-sections, the efficiencies 
and the resolutions and can be assumed to be normalized to 1 without loss of generality but 
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in general do not have the same shape[]. The total asymmetry a*°* receives three different 
contributions. The first comes from the weight functions /, / and / 7^ / corresponds to 
an asymmetry in the experimental apparatus. If matter effects are relevant then there is a 
second experimental asymmetry to be distinguished from the third "intrinsic" asymmetry 
due to CP-violation. Matter effects and / 7^ /, i.e. the CP-violation of the experimental 
setup, give thus rise to a non-vanishing asymmetry a** p even when the mixing matrix U 
is real. It is thus very important to distinguish the experimental asymmetry from the 
asymmetry a£j p due to intrinsic CP-violation which can be written as [0 

ag p ee a« - a- p . (12) 

By definition, ag p must vanish upon setting all CP- violating phases to zero and measuring a 
non-vanishing ag p would be a signal of intrinsic CP-violation in the vacuum mixing matrix 
U. Extraction of afj from data requires of course knowledge of a^ p which is only possible 
if the CP-conserving parameters involved in neutrino oscillation are known. This is most 
likely the case once a measurement of CP-violation in neutrino oscillation becomes feasible. 
The uncertainties of cl* p will however add to the uncertainties on a£ p . 

The dependence of the total asymmetry on the weight functions and on matter effects is in 
general rather complicated. It is therefore not easy to relate the asymmetry ag p of different 
experimental setups or to a reference setup with vanishing matter effects and CP-invariant 
apparatus, i.e. oscillation in vacuum and the same averaging function / + = (/ + f)/2 
for both channels. For the scenarios considered in this paper we find however that matter 
effects are small enough such that they can be considered as a correction to the case without 
matter. 

We omit therefore from now on the suffix "m" and discuss ag p in vacuum analytically. The 
numerically calculated corrections due to matter effects will be presented afterwards. If 
f± = (f ± f)/2 are the average vacuum weight function and asymmetry, respectively, then 
one finds upon neglecting a twice suppressed J d(L/E)Ppp(u ei — > u ej )(L/E)f-(L/E) in the 
denominator of a]f, 



where 



a^ag p + <f\ (13) 



exp = Jd(L/E)P CP (u et -> u ej ){L/E)f.{L/E) 
^ -fd(L/E)P CP (v e ^v ej )(L/E)f + (L/E)> 1 &j 



Jd(L/E)PM^ - v e ML/E)f + {L/E) 



« Jd(L/E)P CP (u ei ^u ej )(L/E)U(L/E)- 1 1 

The quantity a°J p is now the remaining contribution to a*°* coming from the experimental 
apparatus, or more precisely the asymmetry averaged with the CP-conserving part of the 
probability, whereas af? is the CP- violating contribution to a^, given by the CP- asymmetry 



5 The original weight function fo(L,E) is given by the experimental energy and length spectrum and 
depends thus on L and E. Since P is a function of L/E only we can introduce the length averaged weight 



function f(L/E) = / (E 2 /L) f (L,E/L ■ L)dL. 
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averaged with /+. A nice feature of the asymmetry a^ p in vacuum is that it does not depend 
too much on the function / + , namely on the details of the experimental apparatus. Most of 
this dependence cancels in fact in the ratio, especially when the two frequencies giving rise 
to the asymmetry are well separated. |a^ P | < 1 follows from \P&p\ < Pcp- We will study 
in the following first the CP-violating contribution a^ p in vacuum, assuming that a*J p is 
small or under control. Matter corrections will be presented subsequently. 

2.4 Existing experimental data 

Today, there are three different kinds of experiments, which are in favour of neutrino os- 
cillations. First there is the long known solar neutrino-deficit which can be explained by 
an oscillation u e — » v x . In the context of the standard solar model ||, the data, which 
are mainly obtained from Super-Kamiokande, SAGE, Gallium and Chlorine experiments, 
can be fitted by a two neutrino oscillation either in vacuum or matter enhanced (MSW- 
effect). The vacuum solution requires a large mixing with Am 2 m (0.5 — 2.0) • 10~ 9 eV 2 . 
For the MSW-solution we have three different allowed parameter regions [|J where the low- 
Am 2 solution is somewhat disfavoured compared to the solutions with high Am 2 , which 
split into a small mixing (SMA) and large mixing (LMA) case. The SMA solution with 
Am 2 « (0.4 - 1.0) • 10~ 5 eV 2 and the LMA solution with Am 2 » (0.2 - 2.0) • 10" 4 eV 2 (at 
99% C.L.) are obtained from global fits of the rates measured by all solar experiments and 
the day/night measurements of Super-Kamiokande 0, ||. The LMA solution^ will be espe- 
cially interesting in the case of three neutrinos from the point of CP-violating effects, while 
scenarios with more neutrinos have interesting CP-violating effects even without the LMA 
solution. For CP-violating effects in connection with the LMA-solution it is also important 
that the upper bound for Am 2 depends especially on the Chlorine experiment |J. As a 
consequence, excluding or weakening the Chlorine results in the analysis can significantly 
enlarge the possible CP-violating effects. 

The second evidence for neutrino oscillations is given by the atmospheric neutrino data. 
The Super-Kamiokande experiment measures a zenith angle dependent flux of atmospheric 
muon neutrinos which can be explained by an oscillation of the type — > i/ TjS with a 
value of Am 2 « (0.3 - 8) • 10~ 3 eV 2 (at 99% C.L.) and maximal mixing; M. The ratio of 
neutrinos reaching the detector to the number of neutrinos produced in the atmosphere, 
Rft,e = N(i^, e + z^u,e)/No(*V + ^,e), is given by 

Rp = P(vp -> Vp) + RQ X P C p{v e -> Vp) + RQ l alPop{v e -> , (15a) 
R e = P{v e -> u e ) + R Pcp(ve -> Vn) - Ro a oP&p(v e -> Up) , (15b) 

with the particle-antiparticle asymmetry a^' 6 = iVo^^e — v^^)/Nq{v^^ + u^e) in the initial 
neutrino flux and the initial electron-muon asymmetry Rq = Nq{p i1 + z/ At )/Ao(t / e + ^e)- 
Super-Kamiokande measures a zenith angle dependence (and therefore a L/E dependence) 
of R^, but none for R e . 

6 If the last data points of the electron recoil energy spectrum in Super-Kamiokande are reliable, then 
they would tend to favour the SMA solution and to disfavour the LMA solution. This is however still largely 
an open problem and we will delay this issue therefore. 
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A third indication for oscillation is claimed by the LSND experiment ||. It has to be 
remarked, however, that major parts of the LSND allowed parameter region are in contra- 
diction to the KARMEN experiment. This results needs to be checked therefore by new 
experiments (e.g. MiniBooNE). If the LSND result is correct, then it implies oscillations of 
the type — > v e with a lower limit of 10 _1 eV 2 for the squared mass difference and good 
fits around one eV 2 . 

Apart from these positive results, important constraints for the oscillation parameters are 
provided by negative results of disappearance experiments of the type v e — > v e and — > 
v p. H> H3- An important result stems from the Chooz experiment, which allows to put 
an upper limit of about 10 _3 eV 2 on all frequencies contributing to v e disappearance with 
an amplitude larger than 0.2. The Bugey experiment found strong constraints JTTJ] on the 
amplitudes of contributions to v e disappearance with frequencies larger than 3 ■ 10~ 2 eV 2 . 
On the other hand the CDHS and CCFR experiments [|l|, |H| give a limit on the amplitudes 
of contributions to disappearance with frequencies larger than 0.3 eV 2 . There exist also 
absolute neutrino mass limits, like for example from precise measurements of the endpoint 
in the /?-decay spectrum, which lead in principle to further upper bounds on Am 2 . These 
limits are currently however weaker than the ranges quoted above. 

3 CP-violation with three neutrinos 

As already mentioned it is not possible to accommodate the three different experimental 
signals for oscillation in scenarios with only three neutrinos, with three different squared 
mass differences, |Am 2 „-|, i < j, among which only two are independent. One must therefore 
either exclude one evidence for oscillation from the analysis or postulate the existence of 
further neutrinos. Since these options are very different we will study here both possibilities. 

The LSND evidence for oscillation is in a large part of the allowed parameter space in 
contradiction with limits from KARMEN. The LSND evidence for oscillation is therefore 
by far the most controversial one, while the atmospheric and solar ones seem more solid. 
We will therefore consider in this section first the case in which there are only three light 
neutrinos and where the LSND evidence is left out. In the next section we study a four 
neutrino scenario including also the result from LSND. Both studies will include all further 
relevant exclusion limits from experiments with negative results and we will see that CP- 
violation effects can in both cases be quite sizable. 

In the following discussion of the three neutrino scenario we will call v\ and v 2 the mass 
eigenstates which correspond to the smallest |Am 2 |. The ATM and SUN results imply a 
hierarchy between the relevant squared mass differences^ which leads to |Am 2 2 | <C |Am 2 3 | ~ 
|Am 2 3 |. In this limit we can write for eq. (§) 

P C p{v ei - ^) ^ Sij -A\U l3 \ 2 (S lJ - \U j3 \ 2 )sm 2 A 23 -ARe(U jl Ul i U l2 Ul J )sm 2 A 12 , (16a) 

P&p{v £i -> v e .) ~ 8cr i:) Jsin 2 A 23 sin A12 . (16b) 

7 This hierarchy of frequencies can correspond to a situation in which one mass is much larger than the 
other two, ma mi,m2 and in particular, but not necessarily, S> mi 3> mi, to a situation in which 
two masses are degenerate, m\ ~ mi ^ or to a completely degenerate situation mi = m2 ~ m.3. 



8 



The implications of solar and atmospheric neutrino experiments on the parameter space 
are easy to recover due to the constraints on the matrix element U e z between the electron 
neutrino and the third mass eigenstate. First of all a large \U e 3\ 2 (i.e. close to unity) 
is excluded because, due to unitarity, it would prevent solar neutrinos from oscillating. 
Moreover a three neutrino fit of the atmospheric neutrino data |14] gives \U e s\ 2 < 0.08 
at 95% CL. Finally if |Am| 3 | > 2 ■ 10~ 3 eV 2 the results of the Chooz experiment || give 
| C/ e 3 1 2 < 0.05 at 90% CL. These constraints are strong enough to decouple the solar and 
atmospheric neutrino analysis. In fact in the limit of small |?7 e 3| 2 , namely small #13 in the 
parameterization of eq. @, it is easy to recover from eqs. (|Tj]) that the oscillation probability 
for solar neutrinos isQ 

P{y e -»• v e ) ~ 1 - A\U el U e2 \ 2 sin 2 A 12 ~ 1 - sin 2 29 12 sin 2 A 12 , (17) 

so that the solar neutrino plots have to be read in the sin 2 2#i2-| Am 2 2 | plane. On the other 
hand the probabilities involved in atmospheric neutrino experiments are P[y e — * ^e) — 1> 
P(z/ e z/ M ) ~ P(vp -> v e ) ~ , 



P(^ -> ^) ~ 1-4^1(1 - |^ 3 |rsin 2 A 23 ~ 1 - sin 2 2^3 sin 2 A 
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so that the atmospheric neutrino plots have to be read in the sin 2 2^23^1 Am 2 3 | plane. There- 
fore 0, 13 from the experimental data follows 



|Am 23 | = Aibatm = (0.3-8) • 10~ 3 eV 2 atmospheric (ATM) , (19a) 

r (0.2-2.0) • 10" 4 eV 2 large angle MSW (LMA) 
|Am 2 2 | = Amg UN = I (0.4-1.0) • 10- 5 eV 2 small angle MSW (SMA) , (19b) 
{ (0.5-2.0) • 10^ 10 eV 2 vacuum 

as well as sin 2 2#23 ^ 0.8 and sin 2 26*i2 > 0.7 for the large angle MSW solution, all of them 
at 99% CL. For given Am^ we can see from eq. (|16|b) that P&p becomes maximal if Am 2 2 
is in the MSW range for maximal J, i.e. maximal mixing angle. This is the LMA case 
and we can see now why it is especially interesting from the point of CP-violating effects in 
the case of three neutrinos. The point is that this solution is neither affected by the small 
angle suppression of the SMA solution nor by the small Amg UN suppression of the vacuum 
solution. 

The CP- violating asymmetry a^ p , which depends on the parameters |Am 2 2 |, jAm^l, #12, 
#23 ; #13 and 5 in different ways, can now be analyzed. The largest impact on a^ p comes 
from 613 and |Am 2 2 | since 813 can have any value between zero and the limit discussed 
above and |Am 2 2 | is the most important suppression factor. |Am 23 | is also important since 
it controls the relative importance of the sin 2 A 12 terms compared to the sin 2 A 2 3 terms in 
oscillation probabilities. On the other hand 6 I 12 and #23 do not affect a^ p very much since 
they are rather strongly constrained. Finally the (^-dependence of a^ p is of course crucial. 
a^ p is however in leading approximation only linear in sin 5 and it makes therefore sensed 



8 Eq. (0) holds in vacuum, but the argument is essentially unchanged in presence of matter effects. 
9 The quantity \a^ Vc / sin <5| is independent of <5 only in the limit in which the following eq. ( pl| ) holds. 
The lower parts of the contour plots of |a^ e / sin <5| in fig. |l] have in fact a dependence on 5, whereas the 
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to study a^ p / sin 5 instead of a^ p . We will therefore discuss for the rest of this section in a 
quantitative and qualitative way the dependence of |a^ p / sin<5| on 9 13 , |Am 2 2 | and |Am| 3 | 
in the channels — > v^jv^ — > z/ e , — > v T /v^ — > £v, of long-baseline experiments. These 
are in the case of three neutrinos the only experiments where CP-violating effects have a 
chance of being sizable. Short-baseline experiments can not be affected by CP-violating 
effects since the largest squared mass difference is of order 10 _3 eV 2 . Atmospheric neutrino 
experiments are in the case of three neutrinos (and we will see also in the four neutrino 
case) also not much affected by CP- violation, as we will explicitely see in a moment. The 
smallness of matter effects allows us, as mentioned before, to analyze the asymmetries first 
analytically in vacuum. A numerical determination of the matter effect corrections will 
follow at the end of this section. 

3.1 Long-baseline —>■ v e 

Let us consider the — > v e channel and study first the qualitative features of a^ p for the 
large mixing angle solution. For sin 2 2# 12 = sin 2 2# 2 3 = 1 one obtains 

sin 2 26>i3 (sin 2 A 23 ) 
Pop K - v e ) ^ ^ 11 , (20a) 

\P&p{vii — > v e ) \ = I cos 6*i3 sin 26*13 sin 5 (sin 2 A 2 3 sin A i2 ) | . (20b) 

Even though a long-baseline experiment is not sensitive to the suppressed |Am 2 2 | terms 
through the CP-conserving part of the probability (which is quadratic in sinA^), it can 
still be sensitive to |Am 2 2 | through the CP- violating term for two reasons: First P&p is 
only linearly suppressed by sinA 12 . Second Pep is suppressed by sin 2 2#i 3 (which can 
not be large), whereas Pqp contains only sin 26*13, so that experiments able to detect a 
CP-conserving probability which is twice suppressed by sin #13 can see larger asymmetries. 
From eqs. Q20) follows 



sin 2 A 23 sin A12) . 

sin 



sin 2^ (sin 2 A 23 ) 



(21) 



that shows how the (sin A12) suppression is balanced by the 1/ sin 26*13 enhancement of a„ V( . 
The weight function /+ (i.e. experimental details) enters in the asymmetry a CP to a good 
approximation only via (L/E) in (sinAi 2 ) ~ Am 2 2 /4 (L/E). 



in 



23 I 



Fig. [I] shows the contour lines for \a„ V J sin 5\ in the | Amf 2 |-sin 26*13 plane for fixed | A 
(fig. |l|a) and in the |Am 23 | -sin 2 2^13 plane for fixed |Am 2 2 | (fig. |l|b,c,d) using the exact 
formulas for the oscillation probability in eqs. (B) and setting sin 26*12 = sin 2 2023 = 1, which 
as mentioned does not affect the generality of our results. We assumed an experiment with a 
baseline of L = 730km and a energy distribution around E = 6 GeV with <je ~ 3 GeV, giving 

upper parts are almost (^-independent. Fig. [j] assumes | sin <5| = 1. For smaller values of |sin<5| the lower 
parts of the plots in fig. |l| correspond to larger or smaller values of |a^ e / sin <5| according to the sign of 
cos 5 but at the same time la^J gets smaller. When |sin5| = 1 the plots do not either depend on the 
angle 612 corresponding to sin 2 2# 12 = 1 or on the angle #23 corresponding to sin 2 26*23 = 1- 
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(L/E) ~ 100m/ MeV. This looks like the MINOS @ setupQ, but the results depend to a 



good approximation only on {L/E). For other experiments with different values of (L/E) 
one can therefore rescale the asymmetries of fig. [I] with a factor (L/E) /(100m/ MeV). The 
horizontal shadowed regions limit the range of | Am 2 | according to the values given in eq. (|T9"D 



(and |Am 23 | > |Am 2 2 | in fig. |l|d). The MINOS sensitivity taken from ref. (thick solid 
line) and the region excluded by Chooz and the atmospheric neutrino fits (vertically shaded 
regions) are also displayed. The figures show that the parameter space which is accessible by 
the MINOS experiment is not that large. Nevertheless there could be maximally a 30-40% 
effect. 

The general structure of figs. [I] can be easily understood with the help of the approximations 
yielding eq. (|2l|) . Fig. |l]a shows that for maximal CP- violation (| sin<5| = 1) the asymmetry 



can easily reach rather large values. Of course larger values of |Am 2 2 | are preferred. As long 
as sin 2 2#i 3 is not too small (otherwise the sin 2 A23 terms in eq. (||b) become as small as the 
sin 2 A 12 ones and eqs. (^) is no longer valid) also smaller values of sin 2 2# 13 are preferred, 
as expected. Figs. [l]b,c,d show the dependence on |Am| 3 | in greater detail. 

As long as the approximations leading to eq. Q2"T| ) hold, \a^ p Ue / sin5\ oc (sin 2 2^ 13 )" 1 / 2 for 
fixed |Am 2 2 | and the contours are vertical a part from the (small) effects of sin 2 A 2 3 in 
^sin 2 A 23 sin A 12 ) / ^sin 2 A 23 ). On the other hand, for lower values of |Am 23 | the sin 2 A 12 
terms in the CP-conserving part of the probability are not negligible anymore (otherwise 
eq. ( [21~D would give |a^J > 1 for sin 2 2#i 3 low enough). The value of |Am 23 | at which the 
sin 2 Ai 2 terms fold the contour-lines is higher when the sin 2 A 23 term in Pep is smaller, 
namely when sin 2 2#i 3 is smaller, on the left-side of the plot. 

In fig. |l|b a "CP-disfavouring" value of |Am 2 2 | = 0.5 • 10~ 4 has been chosen, while in 
fig. [T|c the "CP-optimistic" value |Am 2 2 | = 2 • 10~ 4 has been used. In fig. ||d a value of 
|Am 2 2 | possible only in some non-standard solar analysis (see below) and therefore lying 
in the shadowed region in fig. |l|a has been considered: |Am 2 2 | = 8 • 10~ 4 . We see from 
figs. |T]b,c that in case of maximal CP-violation the effects in the allowed region covered by 
the MINOS sensitivity range from 2% to 10% (|Am 2 2 | = 0.5 ■ 10" 4 eV 2 ) and from 10% to 
40% (|Am 2 2 | = 2 ■ 10~ 4 eV 2 ), according to the value of |Am 23 | and the amplitude of the 
— > v e oscillation. 

Even in the framework of a standard analysis of solar data, these effects can therefore be large 
enough to spoil the analysis of a possible signal measured by a long-baseline experiment, if 
this analysis does not take into account CP-violation. In figs [l] this can be seen explicitely 
from the shown sensitivity goals of MINOS. Regarding the possibility of measuring such an 
asymmetry, it should be noticed that the asymmetry is larger when the amplitude is smaller, 
so that an enhancement of the CP-violation effect is accompanied by a suppression of the 
statistics. Thus there is no advantage from a statistical point of view. Large asymmetries 
should make it however easier to distinguish intrinsic from experimental asymmetries. 

In some versions of solar neutrino data analysis |Am 2 2 | can lie in the lower part of the 
atmospheric range. This can happen, for example, if one assumes an unknown large sys- 

10 The MINOS experiment can run in different configurations. The energy distribution used here refers to 
the so called "PH2 (medium)" initial neutrino flux, since it gives a sensitivity comparable to the PH2(high) 
flux but a better value of (L/E). 



11 




Figure 1: Contour lines for |a^ e /sin5| in the sin 2 29 13 -|Am 2 2 | plane for |Am| 3 | = 
2 • 10" 3 eV 2 (a) and in the sin 2 20 i3 -| Am 2 3 | plane for |Am 2 2 | = 0.5 • 10" 4 eV 2 ,2 • 
10" 4 eV 2 ,8 • 10" 4 eV 2 (b,c,d) for sin 2 26 l2 = sin 2 29 23 = 1 and (L/E) ~ 100m/ MeV (see 
text). The horizontal shadowed regions limit the range of |Am 2 | according to eqs. ftTQ ) (and 
|Am 23 | > |Am 2 2 | in (d)) whereas the shadowed regions on the right correspond to the com- 
bined Chooz and atmospheric neutrino limits. The thick solid line represents the MINOS 
sensitivity. 



12 



tematic error in the Chlorine experiment ||. This is interesting, since the exclusion of the 
Chlorine data from the analysis makes the remaining solar data also consistent with an 
energy independent reduction of the solar flux as it happens to be above the MSW range 
and for almost maximal oscillation amplitude. As a consequence, for sin 2 2# 13 < 0.25, as 
given by the Super-Kamiokande atmospheric analysis, the large angle solution becomes a 
vertical strip in the | Am 2 2 |-sin 2 26*i 2 plane close to the sin 2 2#i 2 = 1 axis and extending 
to the Chooz limitQ |Am 2 2 | < 10~ 3 eV 2 . This explains why we considered in fig. [TJd the 
possibility |Am 2 2 | = 8- 10~ 4 eV 2 , where the asymmetry is always large for maximal phases^. 

It is interesting to observe that the CP-violating part in the — > v e amplitudes at long- 
baseline experiments can provide information on parameters typically accessible to solar 
neutrino experiments (at least in the present 3-neutrino scenario), while the CP-conserving 
part of the amplitudes are almost insensitive to the same parameters. The detection of a 
CP- violation effect would for example select the large angle MSW solution out of the three 
possible solutions of the solar neutrino problem. The values of |a£%J and sin 2 26 l 1 3, from 
the CP-conserving part, would select a range for |Am 2 2 |. This can be understood with the 
help of fig. [TJa. The measurement of sin 2 2# 13 would select a vertical line in that plot. On 
the other hand, since that figure is "maximal", namely plotted for |sin<5| = 1, values of 
|a^,J in this figure lower than the measured one would be ruled out, allowing thus only 
a limited range for |Am 2 2 |. Since lower values of |Am 2 2 | would give a small CP- violation 
effect, such a measurement would select the upper part of the |Am 2 2 | range provided by the 
solar analysis and larger asymmetries would correspond to stronger lower limits on |Am 2 2 |. 

The parameters of our typical long baseline experiment are not too far away from the MSW 
resonance region. We must therefore include matter effects in our discussion as outlined in 
subsection |2.3| and we will now demonstrate that the corrections to the results obtained so far 
are moderate. Matter effects lead as usual to an inherent sensitivity of the charged current 
interactions of neutrinos to the electron flavour and diagonalization of the Hamiltonian leads 
to new mass eigenstates and shifted masses for the propagation in matter. The transition 
from vacuum to matter can altogether be phrased as two parameter mappings 

(012, fis, 023, S, Am 2 2 , Am 2 3 ) - (6' 12 , 6[ 3 , 9' 23 , 6', (Am 2 2 )', (Am 2 3 )') , (22) 

one for neutrinos and one for antineutrinos. In order to determine the corrections due to 
matter effects we studied for f = f numerically the asymmetry a^f for the following cases: 
a) ar£ = a^°* in vacuum for maximal CP-phase, b) in matter with maximal phase, c) 
a e ^P in matter and d) a^f in matter with maximal phase. These asymmetries are plotted 
in figs. 0a-d for the — > u e channel, where we use the scenario previously assumed for 
fig. |l]c, namely a LBL experiment with (L/E) ~ 100m/ MeV, 5m\ 2 = 2 • 10~ 4 eV 2 and 
sin 2 26* 12 = sin 2 26^3 = 1. Figure |2|b shows the CP-asymmetry measured in the experiment, 

11 Whcn |Amf 2 | approaches the Super-Kamiokande atmospheric range, the |Am 2 2 | effects become im- 
portant in the atmospheric analysis so that the constraint sin 2 26\^ < 0.25, as well as sin 2 26*23 — 1 is not 
reliable anymore. 

12 A non-standard solar analysis leaves open the possibility that |Am 2 2 | falls in the atmospheric range 
while remaining under the Chooz constraint. One can then wonder whether i) |Am 2 2 | could be responsible 
for both the solar (Chlorine excluded) and atmospheric evidences with |Am 23 | above the atmospheric range 
and whether ii) |Am| 3 | could explain the LSND signal. The answers are that, independently of ii), i) is 
strongly disfavoured and in any case ii) is not possible. See appendix A 
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Figure 2: CP- asymmetries for the four cases discussed in the text in analogy to fig. 1: (a) 
a^f = a}^ in vacuum for maximal CP-phase, (b) a^° l in matter with maximal phase, (c) 
a^ p in matter and (d) a^f in matter with maximal phase. 
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plot (c) the matter induced asymmetry and (d) the experimental asymmetry with the matter 
induced asymmetry subtracted. For a few parameter points these asymmetries were already 
calculated in and our figures agree perfectly in those points. Note that matter effects 
depend both on the sign of Am 2 and the sign of the CP-phase 5. In fig. ||we used a positive 
5 and a positive Am 2 , in which case matter and intrinsic effects go into the same direction. 

Altogether we can see that matter effects require a generalization of the asymmetry in order 
to isolate genuine CP-violating effects from matter induced effects. Figures |2| show that the 
corrections are altogether moderate. Ideally one would like to study case d) which is due 
to intrinsic CP-violating effects, but would require to determine somehow a e ^P in matter. 
With good knowledge of the masses and mixing angles this may for example be possible by 
calculating this matter-induced asymmetry theoretically. To do so one has to use fig. 0b 
and subtract (with the correct sign) the theoretically calculated fig. |2|c, which leads to 
larger systematic uncertainties. For fig. |||b the region with large CP-effects shifts to larger 
sin 2 2# 13 and smaller Am| 3 and the maximal CP-asymmetry expected in the sensitivity 
range of MINOS increases from 40% to 50%. This shows that the seperation is in principle 
possible, but it is clearly a difficult task which requires sufficient experimental information. 
Another promising method to seperate the matter-induced and intrinsic CP-asymmetries 
by using envelope patterns of the oscillation is discussed by Arafune, Koike and Sato . 

3.2 Long-baseline — ► v T 

In the — > u T channel the CP- violating probability is the same as in the — *> v e channel 
(a part from the sign) whereas the CP-conserving probability is not suppressed by sin 2 2# 13 
anymore and is therefore larger. Therefore the asymmetry is smaller in this channel. On one 
hand the enhancement of the CP-conserving probability gives better statistics and hence in 
principle the possibility to appreciate a smaller asymmetry. This enhancement would on the 
other hand also enhance the "experimental" contribution a cxp to the asymmetry, making it 
very hard to identify the CP- violating contribution a CP . The —>■ v T channel is therefore 
essentially unsuitable for the detection of CP-violating effects, and one can neglected them 
in the analysis of a possible signal. This statement is confirmed by fig. || where the contour 
lines for \o^ p Vr / sinS\ are plotted in the | Am 2 2 |-sin 2 2#i 3 plane for |Amf 3 | = 2 • 10 _3 eV 2 (a) 
and in the | Am 2 3 |-sin 2 2#i 3 plane for the optimistic case |Am 2 2 | = 2 • 10~ 4 eV 2 (b). The 
figures show that the asymmetries are always smaller than 5%. Matter effects may then (via 
its asymmetric influnece on the mass eigenvalues) even dominate the asymmetry. Therefore 
we do not consider the /i-r channel any further. 

4 Scenarios with four neutrinos 

In order to accommodate the LSND signal in addition to the solar and atmospheric results 
we will study in this section four neutrinos scenarios. We expect that matter effects in 
long baseline experiments are moderate corrections like in the case of three neutrinos. The 
discussion of matter effects will therefore be covered elsewhere. Let us first consider the 
possible mass ordering schemes resulting from the hierarchical values of the squared mass 
differences Am| UN <C Am\ TM <c Am 2 SND . As in the case of three neutrinos we associate the 



15 




Figure 3: Contour lines for \a„ p UT / sm6\ in the |Am 2 2 |-sin 2 2# 13 plane for |Am| 3 | = 
2 • l(T 3 eV 2 (a) and in the |Am| 3 |^in 2 26 13 plane for \Am 2 12 \ = 2 ■ l(T 4 eV 2 (b) for 
sin 2 2^i2 = sin 2 2^23 = 1 and (L/E) ~ 290m/ MeV (see text). The horizontal shadowed re- 
gions limit the range of \Am 2 \ according to eqs. flTQ ) (and |Am| 3 | > |Am 2 2 | in (d)) whereas 
the shadowed regions on the right correspond to the combined Chooz and atmospheric neu- 
trino limit. 



smallest squared mass difference with the first two mass eigenstates, i.e. |Am 2 2 | = Am| UN . 
Then there are only two possibilities for the ordering of the other mass differences relative 
to these first two states: 



A) The intermediate squared mass difference occurs between the 3rd and 4th eigenstates, 
i.e. |Am 34 | = Am^ TM and the larger LSND value defines the splitting between the 
lst/2nd and 3rd/4th mass eigenstate. Altogether this implies |Am 2 2 | = Am| UN <C 
|AmJ 4 | = Am 2 ATM < |Am 2 4 | ~ |Am 2 4 | ~ |Am 2 3 | ~ |Am 2 3 | = Am 2 SND 

B) The intermediate squared mass difference occurs between one of the 1st and 2nd 
eigenstate and one eigenstate out of the 3rd and 4th. Conventionally they are the 2st 
and 3rd mass eigenstate, i.e. |Am 23 | = Am\ TM and the larger LSND value defines the 
splitting between the 4th mass eigenstate and the others, i.e. |Am 34 | = Am 2 SND . We 
have thus in this case |Am 2 2 | = Am| UN <^ |Am 2 3 | ~ |Am 23 | = Am\ TU <^ |Am 2 4 | ~ 
\Am 2 24 \ ~ |Am 2 4 | = Am^ SND 



Note that only Am 2 enters in neutrino oscillation experiments and that this leaves some 
freedom in the ordering and absolute values of masses. Scheme B turns out to be in dis- 



agreement with experimental data fL6 |. We will therefore only consider scheme A in the 



following. This can be understood in a simplified picture where only two neutrino mass 
eigenstates (i.e. their Am 2 ) participate in each oscillation experiment together with the 
information about the involved flavour transitions of each experiment. If one starts in sce- 
nario B with Am| UN as the smallest quadratic mass splitting which involves v e and assumes 
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that Am\ TU (which fixes u^) comes next, then the third and largest AmL SND could not 
be any longer an oscillation between and z/ e , which is a contradiction with the LSND 
experiment. 

Scenarios with four neutrinos involve in general a larger number of parameters in the mixing 
matrix with considerably more complexity in the parameter restrictions. The observed mass 
hierarchies allow however in experiments sensitive to Am 2 > 10~ 3 eV 2 the approximation 
|Am 2 2 | — 0, unless |Am 2 2 | is at the upper border or beyond its standard range, which 
will not be considered here. This approximation simplifies the general task considerably 
and reduces the number of involved parameters. Unlike the three neutrino case, in which 
|Am 2 2 | could quite safely be neglected in the CP-conserving part of the probabilities while 
it was crucial in the CP- violating ones, |Am 2 2 | can be safely neglected here completely. The 
oscillation probabilities are in the limit |Am 2 2 | = given by 

Pcp^ -> f«y) = %(1 - 4|f/ i4 | 2 sin 2 A 24 - 4|f/ i3 | 2 sin 2 A 23 ) 

- 4 Re J 3 e f 1 sin 2 A 34 + 4(Re J 3 e f ' + Re J 4 e f 1 ) sin 2 A 24 ( 23a) 

+ 4(Re J 3 e f 1 + Re J 3 T) sin 2 A 23 , 

P&p{v ei -> u e .) = 81m J^ 6l (sin 2 A 23 sin A 34 cosA 34 + sin 2 A 34 sinA 23 cos A 23 ) , (23b) 

where the second and third line of (fffia) are of special interest for our purposes. 

Eq. (p3|b) shows that it is possible to generate the CP-violating part of the probabilities 
from |Am 23 | and |Am 34 |. CP- violation in long-baseline experiments is therefore no longer 
suppressed by the small Amg UN as in the three neutrino case and we will see that it can 
therefore be large. With four neutrinos one can also wonder, whether CP-violation can be 
important in short-baseline experiments able to measure small transition probabilities. We 
will see that, although the effects are not large, a modest improvement of (L/E) would be 
enough to see sizeable effects in the — > v T channel, if the CP-violation phase is large. 
This is because the relative importance of CP-violation becomes larger for smaller effects, 
unlike what happens in the — > v e case, where the relative importance of CP-violation is 
always small. 

Concerning the possibility of CP-violation effects in atmospheric neutrino experiments, we 
notice that in the four neutrino case they are even more unlikely then in the three neutrino 
one. One may wonder why the four neutrino case does not contain the three neutrino sce- 
nario as a specific limit. This is however the case since there is one additional large frequency 
and also further constraints from experiments sensitive to that frequency. Especially impor- 
tant is here the constraint on u e — *> v e from the Bugey experiment which guarantees in this 
scenario that <-> u e oscillations (and therefore CP-violation, which appears only there) 
do not play a role in atmospheric neutrino oscillations. As a consequence CP-violation is 
negligible in atmospheric neutrino oscillations and we will therefore consider in the following 
only long- and short-baseline —>■ v e and —>■ v T experiments. 

From eqs. fl23|) we can see that the oscillation probabilities between two different flavour 
eigenstates v e . and u e . depend only on the 2x2 sub-sector of the mixing matrix involving 
the zth and jth flavour eigenstates and the 3rd and 4th mass eigenstates. That sub-matrix 
is described by 8 real parameters among which 3 are unphysical phases that can be rotated 
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away, one is a physical phase and 4 are mixing parameters. We choose the 5 physical 
parameters as follows: Let v' e . = U i3 v 3 + U^v^ v' e . = Uj 3 v 3 + Uj^u^ be the projections of 
the flavour eigenstates v e . and v e . on the 3-4 mass eigenspace. Then we define analogous 



to ref. |L7| the quantities c, and Cj as the squared lengths of these projections, i.e. q = 
l^l 2 + |£4i| 2 ? Cj = \Ujs\ 2 + \Uj/t\ 2 . Furthermore 9i G [0, 7r/2] is the orientation of v' e . defined 
via ( j 1 , \Ua\) = ^/ci(cos9i,cos9j) and 0^ G [0, 7r/2] is the relative orientation of u' e . and 



v. 



g , cos (pij = \U* 3 Uj 3 + U* A Uj4\/ ' yjCiCj. Finally we define the relevant CP- violating phase as 



5ij = &rg(U* 3 Uj 3 U* 4 Ui4) = arg(J 34 l ). The amplitudes of the oscillating terms in eqs. ( f23| ) 
can be expressed by these parameters to be 

4 Re J 3 ± l = CiCj sin 2^ (cos 20jj sin 26 i + sin 20jj cos 26 i cos Sij) , (24a) 

4 (Re J\f % + Re J 44 ei ) = ^CiCj sin 9i cos 0^- (cos 0;.,- sin 6i + sin 0y cos 6{ cos dy) , (24b) 

4(Re J^ 6 ' + Re J 3 ^ % ) = 4qCj cos ^ cos 0jj (cos 0^ cos ^ — sin 0^ sin 6*j cos Sy) , (24c) 

4 Im J^ 6 * = qcj sin 26* j sin 20^ sin 5^ . (24d) 

The introduced parameters q, Cj, 6i, <pij, can however not be chosen arbitrary since 
the ij/34 sub-matrix is part of a unitary matrix. The parameters must fulfill the following 
unitary constraints: 

< a < 1 ; < cj < 1 , (25a) 
CiCj cos 2 4>ij < (1 — Cj)(l — Cj) . (25b) 

These conditions can be easily understood noticing that the minor can be embedded in an 
unitary matrix if and only if the two C 2 vectors (£^3, U^) and (Uj 3 , Uj±) can be completed to 
a pair of orthonormal C 4 vectors. Eq. (p5|a) corresponds then to the normalization condition 
and eq. (P5a) to the orthogonality condition of the two C 4 vectors. 

The parameters q, Cj, 6^ 0^, are not only constrained by unitarity, but also by the 
results of the v e — > v e and ^ — ► disappearance, atmospheric neutrinos and LSND 
experiments. In order to make these constraints explicit, let us write the formulae for the 
relevant processes in the approximation where |Am| 4 | <C [Am 2 ^! — jAmf.^: 

P(v e -> p e ) ~ 1 - 4c e (l - c e ) sin 2 A 23 - c 2 sin 2 26 e sin 2 A 34 , (26a) 
P{"» -> - 1 - 4c m(! " c m) sin2 A 23 - c 2 sin 2 20„ sin 2 A 34 , (26b) 
-P(^ -> ^ e ) - 4c e c M cos 2 0^ sin 2 A 23 + 2c e c fl sin 20 M sin 20 Me sin 8^ sin 2 A 23 sin A 34 

(26c) 



— c e c^ sin 29 n (cos 20 Me sin 26 ^ + sin 20 Me cos 29 ^ cos 5 Me ) sin A 34 . 

First we notice that only the sin 2 A23 term is relevant in eq. (|26|c) as far as the LSND 
experiment is concerned. The LSND oscillation probability can thus approximately be 
written as P = A -sin 2 A23, where A = Ac e c^ cos 2 Me . The LSND results can thus be plotted 
conveniently in the A-jAm 2 ^! plane. Moreover, the v e — > u e and — > disappearance 
experiments set upper limits both on the amplitudes of sin 2 A23 and sin 2 A3 4 on the r.h.s 
of eqs. (f26|a,b). We are interested in particular in the limits on 4c e (l — c e ) and 4c M (l — c M ), 
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which correspond to two possible ranges for c e and c M , each close to zero or one. For c e , 
however, only the range around zero is allowed since the other close to one would suppress 
solar neutrino oscillations in an unacceptable way. Therefore we have a limit on c e in the 
form < c e < a° e |16|]. Concerning the range of c M , the sin 2 A 2 3 term in eq. (|26|b) is 
approximately constant in atmospheric neutrino experiments besides being constrained by 
disappearance experiments so that c 2 sin 2 29 ^ must account for the zenith angle dependence 
of the measured flux. This selects out of the two possible ranges for the interval close 
to one, that can be written as 1 — a Q < < 1. The other interval around zero would 



furthermore give a too small amplitude for the sin A 2 3 term in eq. (P6|c) accounting for the 



LSND signal, since it would give A < 4aga°, which would exclude the LSND result (see 
fig. fl6| . Finally, since as said c 2 sin 2 29 ^ ~ sin 2 29^ controls the zenith angle dependence 



we also have sin 2 20 ^ ~ 1. The constraint given by the LSND 



of the atmospheric flux, 

experiment on the parameters c e , c^, cos 2 0^ e , namely 4c e c^ cos 2 = A, has not been used 
in previous analysis. As we will see, it will play an important role in the following. The 
solar result does not give any further constraints for this discussion since we are in the limit 
where Am 2 2 = 0. 
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Figure 4: Constraints on the amplitude A = AceC^ cos 2 Me of the LSND transition proba- 
bility. The shown allowed regions in the LSND parameter space correspond to 90% and 99% 
C.L. The upper part of the plot is not shown because it is excluded by BNL E734- The dashed 
lines represent 4a° (lighter) and 4a° (darker), whereas the dotted line represents 4a® a ^. The 
shadowed region is excluded by the condition $H\). 



In the following subsection we will use eqs. (|23| , |24]) and the parameter constraints discussed 
above to study CP-violation in long- and short-baseline experiments. 
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4.1 Long- and short-baseline — > v e 

Let us first identify the allowed parameter range. The viable scheme A) of squared mass dif- 
ferences implies of course |Am 2 4 | = Ami TM = (0.3-9) ■ 10~ 3 eV 2 and \Am 2 23 \ is of course the 
LSND range, jAm^l = Am 2 SND = (0.04-2.8) eV 2 . Furthermore the following constraints 

< c e < a° e 1 — a° < c M < 1 from v e — > v e , — > experiments (27a) 



2 



'fte 



< (1 — c e )(l — c M ) from unitarity (27b) 



4^ cos 2 Me = A from LSND (27c) 

have to be simultaneously satisfied. It is shown in Appendix B that eqs. (^7j) can be fulfilled 
if and only if 

A < min (4aJ(l - aj), 4a°(l - a°)) ~ min (4a°, 4a°) . (28) 

This gives a restriction on the (Am^l range shown in fig. ^ (the shadowed region is excluded, 
the function min(4a° , 4ag) can easily be recovered from the dashed lines). 

From eqs. ( j23| , [2"4] ) it is easy to see that the CP-asymmetries a^ p Ue do not depend explicitely 
on c e , c M . Eq. (|27|c) introduces however a dependence on c e , c M since cos 2 is a function 
of c e , c M . It is therefore important to know the allowed ranges of c e and c M . Eq. (pT|a) alone 
does not guarantee that it is possible to fulfill eqs. (f2~7|b,c). It turns out that it is possible 
to find c M and Me by solving (|27|b,c) if and only if 



A 1 - y/T^A ( n A . 

< c e < min a°, 1 - — . (29) 



4 " 2 - e - V 4a° 



This is a non-wide range when eq. (gg) is fulfilled. For a given c e in the range of eq. (29) 
the possible values of c M are those fulfilling simultaneously 

A A 

— < c u < 1 ; r and c u >l-a°. (30) 

4c e 4(1 c e ) 1 ; 

In this case A/(A.c e c tl ) < 1 and Me is determined by 

cos 2 <p fie = - , (31) 

where c M = 1 can be used as a good approximation^. 

Finally we have all the necessary ingredients to discuss CP- violation in terms of the quantity 
\cLvvJ sin | which depends mostly only on |Am| 4 |, |Amf 3 |, c e , with the ranges given by 
eq. fll9]a), eq. ( p8|) and eq. ( p9] ) respectively. The dependence on sin 2 26^ is in fact negligible 
since sin 2 26 ^ ~ 1 and we can therefore set sin 2 20 ^ = 1 without affecting our results too 
much. Analogously, we set c M = 1. The dependence of \a„J u J sin<5 Me | on 5^ e is also negligible 
in the limit |Am 2 4 | <C |Am| 3 |. 
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More details on the statements of this paragraph can be found in Appendix B. 
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As a consequence of the Bugey limit, the amplitude of the — > z/ e oscillation in a short- 
baseline experiment is necessarily small. One may wonder if such a small oscillation can 
be accompanied by relatively large CP-violation as it can happen in the three neutrino 
scenario. To see that this is not the case it is enough to study the qualitative features of 



l%V sin 5^1 
one gets 



analytically. With the approximations sin A24 — sin A23 and sin 29, 



and in particular 







4 sin Me cos (fi^ 


(sin 2 A 23 sinA 34 ) 


sin 5^ 




4 cos 2 <^ e (sin 2 A 23 ) 


- cos 20 Me (sin 2 A 34 ) 



(32) 







sin 5 Me 





tan 1 



(sin 2 A 23 sin A 34 ) 



'fie 



(sin 2 A 23 ) 



(33) 



for short-baseline experiments. Even in the "CP-optimistic" case with |Am 2 3 | = 8 • 10~ 3 eV 2 
eq. (|33|) is rather strongly suppressed by sin A 34 in a short-baseline experiment. Measurable 
effects could only show up if an enhancement by a large tan</> Me were possible, which would 
require via eq. (|3l|) a rather small value of A/4c e . The smallest value of A allowed by the 
LSND plot and the maximal value of c e allowed by eq. (B^) do however not allow values of 
A/c e much smaller than one. CP- violation is consequently expected to be negligible which 
is confirmed by fig. f)|a, where the contour lines for \o^ v J sin5| are plotted in the c e -|Am 2 3 | 
plane using the exact formulas eqs. (|23| , |2"4|) for the oscillation probability for an optimistic 
value of |Am 2 4 | = 8 • 10~ 3 eV 2 . The assumed experimental setup is L ~ 0.5km and a 
broad distribution for E around 1 GeV which looks like MiniBooNE. Fig. [|a shows that the 
asymmetry does not exceed 3%, even for maximal phase. 

Long-baseline experiments are of course not sin A 34 suppressed and the sin 2 A 23 dependence 
becomes negligible since it is washed out by averaging over the L/E spectrum. Eq. (|32"D 
becomes thus 







sin 5^ e 





2 cos : 



2j3mgS>^ cos 0^ (sin 
,2 



A 34 ) 



■ COS 



'34, 

» /te (sin 2 A : 



^34) 



(34) 



CP-violation in this case is not suppressed at all. This can be seen directly from fig. |b, 
where contour lines for \a^ Ve / sin5 Me | are plotted in the c e -|Am 2 4 | plane for |Am 23 | = leV 2 
in a "MINOS-like" long-baseline experiment like described in the previous section. The 
unshadowed rectangular window in fig. |5|b represents the values of c e which are allowed by 
the Bugey experiment and the unitarity constraint. The CP-asymmetry in the allowed re- 
gion can reach 60% for a maximal CP-violating phase. This leads to the important question 
whether the allowed region could be reached by a long-baseline experiment despite the strong 
constraint from Bugey. The answer depends on the value of |Am 23 | and 5^ e but in many 
cases is yes. A definitive answer would require a plot of the sensitivity of a long-baseline ex- 
periment (here we consider MINOS again) in the allowed region of the parameter space. The 
published sensitivity plots are however for the parameter space {A, Am 2 } of a simple two 
neutrino oscillation in which the transition probability is given by P = A sin 2 Am 2 L/(4£'). 
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Figure 5: Contour lines for \a„^ v J sin 5^ e \ in a four neutrino scenario plotted in the c e - 
|Am2 3 | parameter space of a short-baseline experiment for |Am| 4 | = 8 • 10~ 3 eV 2 (a) and in 
the c e -|Am| 4 | plane of a long-baseline experiment for jAm 2 ^ = leV 2 (b). Only the allowed 
| Am 2 1 ranges are shown. The shadowed regions are excluded by the constraints 



The — > v e transition probability cannot be reduced to that form in our case being rather 

P(vti -» v e ) - -4/2 + c e (-cos20 Me (sin 2 A 34 ) ± sin 2^ sin 5^ (sin A 34 )) , (35) 

where the sign of the (sin A 34 ) term depends on the sign of sin 29^ and can be reabsorbed 
in the definition of 5^. A is the LSND amplitude given in fig. |] and eq. (|2?|c) has been 
used. The sensitivity in the c e -|Am 34 | plane depends therefore on sin 8^, which controls the 
interference in eq. fl35|), and on A (and in turn on |Am| 3 |) mainly through the A/ 2 term. 
Clearly larger values of A are preferred and it turns out that the sensitivity is larger for 
positive values of sin 5. For example the sensitivity obtained for sin 5=1 and | Am| 3 | = 1 eV 2 
allows to reach most of parameter space. Values of A which are in better agreement the 
KARMEN experiment give less sensitivity and the smallest possible value of A allowed by 
LSND at 99% CL would give a transition probability too small to be measured. Nevertheless 
a long-baseline — > v e experiment along the line discussed has in a four neutrino scenario 
a good chance to observe oscillation and this oscillation would most likely contain a sizable 
or even large CP-violating part. 



4.2 Long- and short-baseline z/„ — > v T 



Unlike what we have seen in the three neutrino case, here the — > v T channel turns out to 
be very interesting. From the previous subsection it is immediately clear that the relevant 
parameters are now |Am| 3 |, |Amf 4 |, M , c^, c T , (p^ T and 5^ T . The parameters |Am 2 3 |, 
|Am| 4 |, On and c M are constraint as before. The unitarity constraints eqs. ( p5|) for the /it/34 
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sub-matrix give now additionally the conditions 



C^C r COS 



< c T < 1 , 

< (l-c M )(l-c T ) , 

while is already known to be in the range 

n A . ■ 1 + VI - A 
max(l - oj, ^) ee cJT < c M < . 

In addition to the transition probabilities already given in eqs. ( p6[) there is now the 
channel 



(36a) 
(36b) 



(37) 



P(v^ — > i/ T ) ~ 4c M c T cos 2 sin 2 A 2 3 + 2c M c T sin 26^ sin 20^ T sin 5 MT sin 2 A 2 3 sin A 34 

— c^c T sin 2# M (cos 20^ T sin 26 ^ + sin 20^ r cos 26^ cos 5 MT ) sin 2 A 34 . (38) 

Let us consider first a short-baseline experiment. Among the oscillating terms only sin 2 A23 
can develop and contribute to this transition probability. The r.h.s. of eq. fl3"8] ) is therefore 
to a very good approximation given by the first term which is furthermore suppressed by the 
unitarity constraints on and c T . This suppression turns out to be much less effective in 
the CP- violating part, whose relative importance grows therefore when the total probability 
gets smaller, very much like in the three neutrino — > v e case. One obtains for the short- 
baseline case 







sin 





tan< 



sin 2 A23 sin A 34 ) 



'ilT 



(sin 2 A 23 > 



(39) 



and we can see that the sin A 34 suppression can be compensated by large values of tan 1 



In fig. |a we show contour lines for \o£ p Ut / sin 5 AtT | in the cos 2 - 



'/XT- 



'/IT 



|Am| 4 | 



plane for |Am 



2 1 

23 I 



1 eV for the short baseline experiment described in the previous subsection. The maximum 
possible sensitivity of such an experiment is within the non-shadowed region. The sensitivity 
in less favourable cases is also shown. The precise way how the sensitivity was here obtained 
and the meaning of the "less favourable cases" is explained in greater detail in the following 
long-baseline case. Contrarily to what happened in the three neutrino z/ M — > v e case, in 
which the necessary enhancement was obtained within the sensitivity of the experiment, 
the necessary enhancement is here outside the reach of an experiment like MiniBooNE. 
Note however that an enhancement by a factor ~ 10 in (L/E) ~ 0.5 of the assumed short- 
baseline experiment would be enough to test CP- violation, which can be seen by multiplying 
the asymmetry in fig. |6|a by (L/E) /O.5. It is therefore not necessary to go a long-baseline 
one with (L/E) ~ 100 to test CP- violation. 

Finally we can also discuss in the four neutrino case CP-violation in a long-baseline ex- 
periment like that considered in the three neutrino case. The corresponding contour lines 
for \cLv Ut / sin 8^ are shown in fig. |]b as before in the cos 2 Mr -| Am| 4 | plane for fixed 
|Am 2 3 | = leV 2 . The precise value of lAmSy is anyway irrelevant since sin A 23 is averaged 
to 1/2 in this case. The asymmetry is in this case also completely independent of c M , c T and 



23 



1 0.2 0.04^c r /(c T ) n 



I 10- 3 



: F// 








— 10% : 






— 5% ; 






_ 2.5%; 







10" 4 10" 3 10" 2 10" 1 10° 



cos <p 



> 

1) 



c T /(c T ) max -> 0.05 0.2 



10" 




JUT 



10" 4 10" 3 10" 2 10" 1 10° 
cos 2 0„ T 



Figure 6: Contour lines for \a„ Vt / sin in the cos 2 0^ r -|Am| 4 | plane for |Am| 3 | = 1 eV 
(a) for the short-baseline experiment and (b) for the long-baseline experiment discussed in 
the text. The accessible parameter space of such an experiment depends on c T . The non- 
shadowed region shows the optimal case of c T = c™ ax which covers the whole plot in (b). The 
dashed lines show the sensitivity in less favourable cases with c T < c™ ax (see labels above 
figure and text). Everything to the right of the dashed lines in (a) can be explored by the 
discussed experiment whereas in (b ) the whole area from the upper-left corner to the dashed 
lines lies within the sensitivity of the experiment. 



the dependence on 9^ is negligible as before. What matters is again the region of parameter 
space accessible to the experiment, which depends on c T in the way given by eq. The 
range of c T , which is determined by eqs. (^) in terms of cos 2 ^^, is therefore crucial and 
one finds 

1 — c mm 

< c T < — £-= . (40) 

" - l-c- in sin 2 MT 

The sensitivities of fig. |6|b depend therefore on the chosen value of c T and become maximal 

1 — C m * n 

for the largest possible value c™ ax = - — min M 3 , — . For this maximal value c T = c™ ax one finds 

1 sin (pfir 

that the experimental sensitivity covers the complete coordinate range of fig. |]b and only 
rather small values of c T make some of the relevant parts of the parameter space unaccessible. 
We show in fig. |B|b the sensitivity corresponding to c T = 0.2c™ ax and c T = 0.05 c™ ax . The 
region to the right of the 20% (5%) sensitivity line are excluded if c T is less than 20% 
(5%) of its maximum value in that region. One can see that c T = 0.2 c™ ax still allows to 
reach the interesting region whereas c T = 0.05 c™ ax excludes it. Here, as in the — ► v e 
case, the sensitivity depends on the sign of sin 29 ^ and on the CP- violating phases, that 
control the interference between the CP-violating and CP-conserving parts of the transition 
probabilities. The lines plotted correspond to sin26' At = 1, sin<5 MT = 1, that give a better 
sensitivity. The results of this subsection show that — ► u T is a good place to look for CP- 
violation in four neutrino scenarios even with intermediate-baseline experiments. Moreover, 



24 



if there are four neutrinos and if a signal is observed, CP-violation should be taken into 
account in an analysis. 

5 Discussion and Conclusions 

We studied in this paper CP-violation in neutrino oscillation and discussed the sensitivity 
for such effects of current and future experiments. As a measure of CP- violation, we consid- 
ered asymmetries between CP-conjugated transition probabilities. These asymmetries are 
a measure of the relative strength of CP-violating effects. They are therefore very useful 
in studying how much CP-violation would affect the measurement of oscillation in one of 
the two CP-conjugated channels. The contribution from CP-violation has to be separated 
however from the experimental asymmetries between the two channels, which we assumed 
to be small or under control. 

The fact that it is not possible to accommodate the claimed three independent oscillation 
signals in scenarios with three neutrinos led us to a twofold strategy. The first scenario was 
to leave out one evidence for oscillation and to analyze the case of three neutrinos. Since the 
LSND evidence is almost in contradiction with KARMEN it is considered most controversial 
and we omitted it therefore in the three neutrino case. The second case which was studied 
includes LSND in a four neutrino scenario. In both cases all further existing exclusion limits 
were taken into account. The two scenarios lead to quite different results with different sizes 
of CP-violating effects. In all cases we present results for maximal CP-phase, which can be 
easily rescaled to an arbitrary value. 

For three neutrinos CP-violating effects are drastically suppressed by small angles or by 
extremely small Am 2 in the case of the small mixing angle MSW-solution and for the 
vacuum solution. The large mixing angle MSW-solution allows however sizable CP-violating 
effects in long-baseline experiments, while there is no effect in short-baseline experiments. 
In the /j, — > e channel we found in long-baseline experiments for maximal CP-phase 5 effects 
up to 40%, while we found only small effects in the fi — > r channel. The observation of 
CP-violation would also allow to distinguish between the different solar solutions and can 
even further restrict the parameter space for |Am 2 2 |. Including matter effects in these long 
baseline experiments we found for the considered setup moderate corrections compared to 
the case without matter. Depending on the sign of the squared mass difference, the total 
intrinsic CP-asymmetry in /i — + e oscillation can be enhanced up to 50%. In // — > r 
oscillation matter effects dominate the asymmetry since the intrinsic CP-violation is very 
small and therefore anyway uninteresting. 

In order to include the LSND result we studied also the case of four neutrinos where many 
more parameters exist in principle. The solar Am 2 -value is in this case unimportant since 
CP-violating effects can be generated by the larger Am 2 responsible for the atmospheric 
neutrino oscillations and for the LSND measurement. For experiments which are only 
sensitive to Am 2 > 10 -3 eV 2 , we could make the approximation Am 2 2 = which reduces 
the number of relevant parameters drastically and allows a study of the available parameter 
space. The CP-violating effects are now in the case of four neutrinos potentially larger 
and we considered therefore also short-baseline experiments. Altogether we find in the 
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four neutrino scenario for maximal CP-phase the following effects: In short-baseline fx — > e 
experiments less than 2 %, in short-baseline /i — > r up to 10 % and in long-baseline /i — > e 
as well as /i — > r experiments up to 60 %. The effects are not very big in the considered 
MiniBooNE-like short-baseline setup, but we want to point out, that a modest improvement 
of (L/E) ~ 0.5 by a factor 10 would be enough to test CP-violation in the z/ M — > z/ T channel. 

We did not consider cases with more than four neutrino mass eigenstates. It should however 
be clear from the current analysis that large CP- violating effects could easily be involved in 
that case in current and/or future neutrino oscillation experiments. 

In summary we found that CP-violating effects can be surprisingly large in some future 
neutrino oscillation experiments and such effects should therefore be included in the analysis. 
Besides the obvious general interest for a determination of CP-violation connected to the 
theoretical questions on physics beyond the Standard Model and the potential role which 
CP-violation could play in lepto- and baryogenesis, there are further reasons why CP- 
violation should be taken into account. The main point is that the omission of CP- violation 
can spoil a two or three neutrino analysis that does not take it into account. Moreover, 
if an asymmetry between CP-conjugated transitions were measured and the presence of 
light sterile neutrinos would be excluded, it would discriminate between the different solar 
solutions and set lower bounds on the solar Am 2 . 
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Appendix A 



A non-standard solar analysis allows in principle the possibility that |Am 2 2 | falls in the 
atmospheric range, but still consistent with the Chooz constraint. For this case one may 
wonder whether 

i) |Am 2 2 | could be responsible for both the solar evidences (with Chlorine excluded) 
and the atmospheric results with |Am 23 | above the atmospheric range (to account for 
cosmological requirements on neutrino masses) 

and whether 

ii) |Amf 3 | could explain the LSND signal. 

The answers to these two questions are that, independently of ii), i) is strongly disfavoured 
and that ii) is not possible, i) is disfavoured independently of the possibility of explain- 
ing LSND for two reasons. First of all the probability that |Am 2 2 | is in the part of the 
atmospheric range which is not excluded by the Chooz constraint is small. This, together 
with the solar constraints which exclude large |£7 e3 | 2 gives |£/e 3 | 2 < 0.025. The survival 
probability for solar neutrino state v<i is therefore P{y e — > v^) = 1 — sin 2 2#i2 sin 2 A12 with 
sin 2 29 12 — 1 which has to explain the solar data (with Chlorine excluded) . This is on the 
other hand in this scenario also the survival probability for atmospheric v e appearing in 
eq. (|T5|b). The atmospheric fit would in this case therefore be very bad, since the the zenith 



angle variation of R e in eqs. (|15| ) would be larger than for R^. 

LSND cannot be explained in any case, even if the two points above are ignored. The LSND 
oscillation probability is for this scenario 

P&v -> *e) = 4|t/ e3 ^3| 2 sin 2 A 23 . (41) 

The mixing |t/ e3 | 2 is known to be small with the upper limit given by Bugey. The mixing 
\U^\ 2 is also small because 4|[/ M i[/ Al2 | 2 < (1— \ U^\ 2 ) 2 has to be maximal in order to explain 
the oscillation of in eq. (|T5)a,) and because of the experimental limit given by CDHS 
and CCFR. Combining these two limits into a single limit for 4|[/ e3 £7 Al3 | 2 one finds that the 
probability in eq. fl4lf) is too small to explain the LSND signal 17j . 



Appendix B 

In this appendix we prove the statements concerning the range of c e , c M , Me satisfying 
eqs. (p7|) . First we prove that eq. (|28| ) is sufficient for the existence of a solution to these 
equations. Therefore it is enough to check that c e = c mm , c M = c max , </> Me = are a solution of 
eqs. (|27j). Then we prove that eq. (^) is also a necessary condition. From eqs. (f2~7|b,c) follows 



4c e c M > A and 4(1 - c e )(l - c M ) > A and therefore c e > A/(Ac ll ), c M < 1 - A/(4(l - c e )). 
In particular we have 4c 2 — 4c e + A < and 4c 2 — 4c^ + A < so that, together with the 
previous relations, we obtain 

^ min 1 - y/l-A A max 1 + y/l-A . A 

c e > c mm = ~ — ; c u < c max = ~ 1 . 

2 4 ' M ~ 2 4 
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But in order to satisfy eq. (|2T|a) one must have c mm < a° e and c max > 1 — a? , namely 

A < 4a°(l - a°) and A < 4a°(l - a°) . 

Let us now prove eq. (f29|). To begin, let us consider a given value of c e . Eqs. (|2T|b,c) give 

where a solution exists when eq. (|28| ) holds. Since 1 — a° < c M < 1, 1 — A/(4(l — c e )) must 
be equal to or larger than 1 — a° in order to obtain a finite range for c M . One finds thus 

A 

which, together with the previous bounds, gives eq. (|29|) . To see that eq. (|29l ) is not 



an empty interval we have to check that c™ in < 1 — v4/(4a°). This is a consequence of 
A — (1 — c™ in )(l — c™ ax ) and c™ ax > 1 — a°. The range of c M for a given value of c e , eq. (^0|), 
then follows from what we saw. 
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